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PREFACE 

This  report  is  the  final  technical  report  for  the  Atmospheric 
Phase  Error  Measurements  task  performed  under  Contract  F33657-77-C.-0229. 
The  work  under  this  contract  was  conducted  in  the  Radar  and  Optics 
Division  of  the  Environmental  Research  Institute  of  Michigan  (ERIM) , 

P.0.  Box  8618,  Ann  Arbor,  Michigan.  This  effort  was  sponsored  by  the 
Air  Force  Systems  Command,  Aeronautical  Systems  Division,  Wright- 
Patterson  AFB,  Ohio.  Technical  Monitors  were  E.  Dcardurff,  Jr.,  and 
L.  Faulkner;  the  Contracting  Officer  was  H.R.  Minteer,  ASD/AERK.  The 
Program  Manager  at  ERIM  was  G.F.  Adams. 

The  basic  analysis  was  carried  out  by  R.  Crane  with  data  calcula- 
tions provided  by  H.  Kllmach,  L.  Wiloek,  and  D.  Franczak.  R.  Hamilton 
and  T.  Van  Sickle  assisted  in  documentation.  The  final  report  was 
written  by  G.  Adams  with  major  contributions  from  R.  Crane  and  J. 
Auterman.  Results  of  the  analysis  were  reviewed  by  J.  Auterman, 

R.  Bayma,  R.  Heimiller,  and  L.  Porcello. 

This  analysis  was  planned  jointly  with  M.C.  Thompson  of  the  U.S. 
Department  of  Commerce,  Office  of  Telecommunications.  The  experimental 
data  was  collected  and  provided  to  ERIM  by  the  Department  of  Commerce. 
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INTRODUCTION  AND  SUMMARY 


The  imaging  performance  of  synthetic  aperture  radar  (SAR)  systems 
in  use  today  does  not  appear  to  be  severely  limited  by  atmospheric  pro- 
pagation anomalies.  Atmospheric  limitation  effects  could  become 
important  for  SAR  systems  operating  either  at  longer  ranges  or  with 
finer  resolutions.  Such  systems  would  encounter  an  increase  in  the  amount 
of  atmosphere  through  which  radar  signals  must  propagate  and  then  be 
coherently  piocessed.  This  report  describes  the  results  of  a study  of 
these  atmospheric  effects.  The  study  was  performed  at  the  Environmental 
Research  Institute  of  Michigan  (ERIM)  and  was  carried  out  under 
* the  sponsorship  of  the  U.  S.  Air  Force  Systems  Command,  Aeronautical 

Systems  Division.  It  was  based  upon  experimental  atmospheric  measure- 
« ments  collected  by  the  U.S.  Department  of  Commerce,  Office  of  Telecom- 

muni  cat  ions . 

The  study  which  was  carried  out  is  summarized  in  this  first  section 
of  the  report.  The  theory  of  SAR  systems  is  briefly  reviewed  along  with 
the  nature  of  atmospheric  effects,  then  the  procedures  followed  during 
this  study  are  described,  and  finally  the  results  are  summarized.  The 
d.ita  analvsis  which  was  performed  is  described  in  detail  in  the  second 
section.  Conclusions  and  recommendations  are  included  in  the  final 
section.  Specific  analytical  details  are  included  in  a series  of  appen- 
dices. 

1.1  BACKGROUND  i 

A synthetic  aperture  radar  (SAR)  is  a pulsed  coherent  airborne 
sidelooking  radar  system  which  is  used  to  generate  f ine-resolut ion  radar 
images  of  the  ground.  A SAR  system  images  ground  features  in  a two- 
dimensional  coordinate  system  defined  by  the  nominal  slant  range  from 
the  aircraft  flight  path  and  the  distance  along  the  flight  path.  The 
geometry  of  a SAR  system  is  illustrated  in  Figure  1.  Resolution  along 
the  slant  range  coordinate  is  achieved  through  the  use  of  radar  range 
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pulse  techniques.  Generally  pulse  compression  techniques  are  applied 
using  coded  pulse  waveforms. 

Resolution  in  the  along-track  coordinate  is  achieved  by  coherently 
combining  radar  pulses  received  over  an  along-track  distance  called  the 
synthetic  array  length.  An  antenna  equal  to  this  length  is  synthesized 
and  the  along-track  resolution  is  determined  by  the  synthetic  array 
beamvidth . If  data  from  a constant  real  antenna  beamwidth  is  processed 
for  all  ranges  then  the  synthetic  array  length  increases  linearly  with 
range,  the  resultant  synthetic  array  beamwidth  decreases  with  range, 
and  the  theoretical  synthetic  along-track  resolution  is  independent  of 
range . 

Distortion  of  the  radar  signal  phase  over  the  length  of  the  synthe- 
tic aperture  will  degrade  the  SAR  along-track  impulse  response  performance. 
Therefore  sources  of  sucli  phase  distortions  and  their  effect  upon  SAR 
imaging  performance  is  an  important  consideration  in  the  design  ol  such 
systems.  These  phase  error  sources  can  be  divided  into  two  general 
categories;  those  whose  effect  is  somewhat  controllable  by  engineering 
design  and  those  over  which  the  system  designer  has  no  control.  Exam- 
ples of  the  former  include  system  component  phase  instabilities  and 
uncorrected  aircraft  motion  effects.  Phase  distortions  which  are  caused 
bv  atmospheric  propagation  anomalies  fall  into  the  latter  category. 

These  types  of  phase  errors  cannot,  in  general,  be  corrected  except 
through  the  use  of  sophisticated  adaptive  processing  techniques.  Opti- 
mum performance  in  the  presence  of  propagation  phase  errors  is  an 
important  limitation  of  SAR  performance.  It  is  essential  that  these 
limiting  levels  of  performance  be  understood  and  their  effects  be  in- 
cluded in  system  design  error  budgets.  Thus,  this  study  was  performed 
to  analyze  the  effects  of  such  atmospheric  anomaly  effects. 

1.2  REVIEW  OF  ANALYSIS 

The  analysis  carried  out  during  this  study  is  briefly  outlined  in 
this  section.  The  direct  measurement  of  SAR  signal  phase  errors  caused 
by  atmospheric  anomalies  is  difficult  to  carry  out  in  practice.  The 
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problem  is  one  oi  separating  the  etfects  of  other  phase  error  sources 
from  the  atmospheric  propagation  effects.  Such  phase  error  measurements 
must  he  made  from  an  aircraft  in  order  for  them  to  be  directly  relatnble 
to  a SAK  system.  An  aircraft  cannot  travel  a perfectly  straight  path  at 
a constant  ground  speed.  Yaw,  pitch,  and  roll  variations  as  well  as 
lineai  accelerations  must  he  accurately  determined  and  their  phase  error 
effects  accounted  for  by  the  use  of  a motion  compensation  system.  In 
order  to  measure  the  atmospheric  phase  errors,  the  accuracy  of  the 
motion  compensation  system  must  he  such  that  uncompensated  motion  phase 
errors  are  much  less  than  the  atmospheric  induced  errors.  This  is  very 
difficult,  if  not  impossible,  to  achieve  even  in  a carefully  designed 
experiment.  Therefore,  an  indirect  method  of  determining  atmospheric 
phase  error  effects  upon  SAK  performance  has  been  developed  and  was  used 
during  this  study. 


The  indirect  method  used  is  a relatively  simple  method  based  upon 
the  theory  of  radio  propagation  through  a turbulent  medium  developed 
previously  by  Tatarski  [i|  and  others.  This  theory  was  developed  for 
the  purpose  of  explaining  variations  in  both  the  phase  and  amplitude  of 
radio  waves  due  to  atmospheric  index  of  refraction  variations. 

Thompson,  et  al  [2]  have  shown  that  there  can  be  times  throughout  the 
day  during  which  extremely  large  radio  signal  propagation  amplitude 
variations  can  occur.  At  these  times  the  turbulence  theory  cannot  he 
directly  applied.  However,  during  the  majority  of  time  such  large  vari- 
ations are  absent  and  the  theory  can  be  used  with  a high  degree  ot 
confidence.  In  this  study,  the  Tatarski  theory  is  applied  to  the  analy- 
sis of  SAR  atmospherically  induced  phase  errors.  Amplitude  signal 
variations  are  not  considered  since  the  phase  error  effects  generally 
predominate . 
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Tatarski,  V.I.,  Wave  Propagation  in  a Turbulent  Medium,  hover 
Publications,  Inc.,  New  York,  19M. 


Thompson,  M.C.,  Jr.,  H.B.  Janes,  L.E.  Wood,  and  0.  Smith,  Phase 
and  Amplitude  Scintillations  at  9lb  GHz  on  an  Elevated  Path,  1ITF. 
Trans.,  AP,  November  1975. 
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The  Tatarski  model  does  not  directly  apply  to  the  SAR  problem. 
I’oreello  [1]  earlier  used  the  Tatarski  model  along  with  measurements 
by  Ihompson  and  Janes  [4]  to  determine  the  statistics  of  the  random 
atmospheric  induced  phase  perturbations  in  SAR  signals.  This  earlier 
analysis  by  Porcello  has  been  extended  in  this  report  to  include  the 
simple  first-order  effects  of  a curved,  rather  than  a flat,  earth. 

Brown  and  Riordan  [5]  then  determined  the  average  impulse  (point  target) 
response  of  a SAR  when  phase  errors  caused  by  the  atmosphere  were  pre- 
sent in  the  radar  data. 

Brown  and  Riordan  characterized  the  impulse  response  by  the  use  of 
the  analytical  equivalent  rectangle  measure  of  resolution.  This  measure 
corresponds  directly  to  the  reduction  of  the  mainlobe  response  peak 
amplitude.  It  does  not  differentiate  between  a broadened  mainlobe  and 
increased  side  lobes.  During  this  study  expected  average  impulse  re- 
sponse curves  were  calculated  and  plotted  for  various  levels  of  atmos- 
pheric turbulence.  Then  several  various  measures  of  impulse  response 
quality  were  calculated  for  the  expected  impulse  response  curves.  These 
include  measures  of  mainlobe  width  and  sidelobe  levels. 

A portion  of  the  average  impulse  response  degradation  caused  by 
atmospheric  turbulence  can  be  characterized  as  a simple  quadratic  de- 
focus ing  of  the  radar  image.  The  magnitude  of  this  effect  was  computed 
and  separated  from  the  remaining  higher  order  effects.  Also,  the  linear 
portion  of  the  phase  error  which  causes  a simple  image  translation  was 
computed  and  separated  from  the  phase  error  term. 

The  numerical  results  of  this  study  are  based  upon  atmospheric 


3.  Porcello,  L.J.,  Turbulence-Induced  Phase  Errors  in  Synthetic 
Aperture  Radars,  IEEE  Trans.,  AES,  September  1970. 

4.  Thompson,  M.C.,  Jr.  and  H.B.  Janes,  Measurements  of  Phase-Front 
Distortion  on  an  Elevated  Line-of- Sight  Path,  IEEE  Trans.,  AES, 
September  1970. 

5.  Brown,  W.M.  and  J.F.  Riordan,  Resolution  l.imits  with  Propagation 
Phase  Errors,  IEEE  Trans.,  AES,  September  1970. 


index  ot  retract  inn  prof  11  on  measured  tiy  Thompson,  et  al  |(>|.  The  data 
collected  were  used  to  calculate  values  of  the  Tatarski  model  statisti- 
cal propagation  parameters.  These  values  were  then  combined  with 
various  postulated  radar  cont Igural ions  (ranges,  altitudes,  and  synthe- 
tic aperture  lengths)  to  determine  numerical  estimates  ot  SAK  imaging 
performance  as  limited  by  atmospheric  effects  during  the  days  ot  the 
measurements  at  the  measurement  locations.  The  data  used  were  all 
collected  during  three  brief  series  of  test  flights,  one  over  a 
Florida  coastal  region  and  two  over  a portion  of  Colorado.  These  three 
sets  ot  data  exhibit  different  results  as  would  he  expected.  The  re- 
sults ot  the  calculations  from  these  measurements  are  Included  In  this 
report  . 

1.3  SUMMARY  OF  RESULTS 

The  results  ot  this  study  are  brlelly  summarized  In  this  section. 

A limited  set  ot  experimental  atmospheric  index  ot  retract  ion  data  which 
was  collected  by  Dr.  M.C.  Thompson  ot  the  D.S.  Deparment  ot  Commerce 
was  analyzed  during  this  study.  The  effects  ot  atmosphere  upon  SAK 
imaging  performance  was  determined  bv  extending  previous  analyses  ot 
turbulence  theory,  statistical  theory,  and  SAK  theory.  These  results 
were  then  applied  to  the  data  collected.  It  should  be  noted  here  that 
the  collected  data  was  very  limited  in  quantity.  The  data  was  gathered 
in  limited  geographical  areas,  over  a limited  period  ol  time,  and  on  a 
limited  number  ol  days.  The  capability  ol  drawing  ma|or  conclusions  01 
extrapolating  to  other  conditions  is  therefore  limited. 

Based  on  the  data  collected  and  analyzed  it  was  shown  that  cleat 
air  turbulence  el  feds  would  not  lie  expected  to  seriously  degrade  the 
performance  ol  t he  currently  operational  AN/UI’D-A  and  AN/lll’D-b  SAK 


It.  Thompson,  M.C.,  Jr.,  K.K.  Marler,  and  K.C.  Allen,  Measurement  ot 
Refractive  index  Parameters  for  Calculating  Limits  of  SAR  Image 
Resolution,  OT  Technical  Memorandum  7 7-2  It,  U.S.  Department  ol  Cont- 
meree/Offl ee  of  Telecommunications,  Boulder,  Colorado,  April  l‘*77. 
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systems.  Further  It  was  shown  that  the  performance  of  SAR  systems 
operation  at  significantly  longer  ranges  or  finer  resolutions  may  not 
he  seriously  Limited  hy  atmospheric  effects.  In  fact,  from  the  col- 
lected data  only  moderately  degraded  image  quality  was  calculated  most 
ot  the  time  tor  SAR  systems  operating  at  ranges  as  great  as  100  nmi  and 
Impulse  response  widths  as  fine  as  5 feet.  On  rare  occasions  clear  air 
turbulence  was  shown  to  severely  affect  the  image  quality  through  phase 
error  effects.  It  is  likely  that  more  often  weather  attenuation  and 
baekscatter  effects  as  well  as  motion  compensation  effects  would  have 
seriously  limited  performance. 

The  data  analysis  also  showed  that  there  was  no  discernible  rela- 
tionship between  ground  weather  measurements  and  the  expected  SAR 
performance  during  the  time  when  the  data  was  collected. 

The  data  was  collected  during  a total  of  17  flights  in  two  areas, 
Colorado  and  Florida,  during  the  fall  ol  197b.  Although  the  collected 
data  has  been  extensively  analyzed,  there  is  insufficient  data  to  begin 
Lo  predict  the  limitations  as  a function  of  geographic  location,  season 
time  of  day,  or  weather  condition.  There  is  insufficient  data  to  pre- 
dict expected  performance  variance,  only  the  average  performance  was 
predicted  for  individual  flights.  Also,  reliable  measurements  were  not 
obtained  in  heavy  rain  or  in  heavy  cloud  cover  and  these  cases  are  not 
well  analyzed  to  date.  It  should  be  noted  that  measurements  were 
gathered  in  an  area  expected  to  yield  relatively  severe  effects,  loculi 
adjacent  to  heavy  cloud  cover  and  rain,  and  in  light  cloud  cover  and 
moderate  rainstorms.  It  is  recommended  that  further  data  collection 
and  analysis  be  considered  to  address  these  limitations. 
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ANALYSIS 


The  goal  of  this  analysis  Is  to  determine  the  extent  of  adverse 
effects  on  the  performance  of  a synthetic  aperture  radar  (SAK)  due  to 
the  random  phase  distortions  which  are  Introduced  bv  atmospheric  propa- 
gation anomalies.  The  phase  distortion  is  assumed  to  he  caused  hy 
turbulence  in  the  propagation  medium.  The  analysis  which  was  carried 
out  is  described  here  in  two  separate  parts.  The  tlrst  part.  Section 
2.1,  describes  the  theoretical  model  developed  for  the  analysis  which 
is  based  upon  turbulence  theory  by  Tatarskl  and  others.  The  second 
part,  Section  2.2,  describes  the  data  collection  carried  out  by  Dr. 
Thompson  of  the  Department  of  Commerce,  the  analysis  of  the  collected 
data,  and  the  results  of  that  analysis. 

2.1  ANALYTICAL  MODEL 

The  analytical  model  is  developed  in  a series  of  three  steps.  The 
first  step  of  this  development  is  to  determine  the  form  of  the  varia- 
tions in  tlie  index  of  refraction  of  the  atmosphere.  The  atmosphere  is 
characterized  by  a parameter  called  the  atmospheric  index  of  refraction 
structure  function.  This  function  describes  statistically  how  the  index 
of  refraction  varies  along  different  paths  through  the  atmosphere.  The 
second  step  is  to  calculate  the  effect  of  the  atmosphere  upon  the  phase 
of  SAR  signals  across  the  synthetic  aperture.  This  is  achieved  by  In- 
tegrating the  expected  index  of  refraction  variations  along  radar  rays 
between  different  portions  of  the  synthetic  aperture  and  the  ground  tar- 
get being  imaged.  A SAR  phase  error  structure  function  results  from 
this  step.  This  structure  function  characterizes  the  expected  phase 
variations  across  the  synthetic  aperture.  The  third  step  is  then  to 
examine  the  effect  of  this  expected  phase  error  upon  SAR  image  quality. 
The  expected  SAR  Impulse  response  results  from  this  step.  This  impulse 
response  thus  characterizes  the  average  SAR  performance  as  a function  of 
atmospheric  turbulence  levels  and  system  geometry. 


j . 


8 


A briet  description  ot  t ho  details  and  the  results  ot  each  ot  these 
.malvtlcal  steps  Is  Included  below.  Further  details  are  included  In  the 


appendices  tor  those  portions  ot  the  analysis  which  are  not  general Iv 
available  in  the  referenced  literature. 

.1.1.1  ATNOSrilF.RU'  TURBULKNl’.K  MODKl, 

Tatarski  |ll  has  analyzed  the  problem  of  radio  wave  phase  errors 
caused  bv  propagation  t hrough  a turbulent  medium.  His  analysis  is 
based  upon  a term  called  t he  atmospheric  Index  of  refraction  structure 
function.  This  term  is  used  to  characterize  the  atmosphere  and  Is  de- 
t Ined  as 

n(r)-K|n(x  + r)-ti(x>r  ( 1 ) 

n 

The  parameter  nix')  is  the  Index  ot  refraction  of  the  atmosphere  mea- 
sured along  a straight  path  x.  The  operator  F|  ) denotes  the  expected 
value.  Thus,  this  structure  function  is  the  expected  value  ol  the 
square  ot  the  variation  of  the  Index  of  refraction  over  a distance  r. 

i'atarskl  uses  physical  arguments  to  show  that  for  microwaves  this 
structure  function  is  approximately 

2 2/  1 

1)  (r;  h)  - i'.  ihl  r (2) 

n n 

where  v <<  r <v  1. 

O o 

where  r is  measured  a 1 one,  a horizontal  axis.  ('  (hi  is  the  retractive 

n 

index  constant  along  a radiation  path  and  is  generally  dependent  only 

upon  height  h.  t'he  terms  C and  I.  arc  called  the  inner  and  outer 

o o 

scales  ot  the  turbulence.  The  relationship  shown  for  r and  C , 1.  is 

o o 

generally  assumed  to  be  satisfied  for  the  SAK  cases  which  are  analyzed 
herein.  Figure  illustrates  the  model  geometry. 


2.1.2  SYNTHETIC  APERTURE  PHASE  MODEL 


■KBS*”  l-.-iii’. 


The  structure  function  for  the  phase  distribution  across  a synthe- 
tic aperture  is  defined  as 


D <r) 

41 


E|4>(x  + r)  - 0(x)  1 


2 


f 3) 


where  x is  measured  along  the  aperture, 
by  applying  the  relationship 


The  value  of  D, (r)  can  be  found 

4' 


4>(K) 


R 


n(«)  dC 


(4) 


where  k » 2n/X 

R * path  length. 

The  factor  2 must  be  Included  since  the  radar  signal  travels  along  the 
antenna  to  target  propagation  path  twice,  once  to  illuminate  the  target 
and  once  as  reflected  by  the  target. 

The  procedure  followed  herein  is  to  determine  the  effective  phase 

distribution  across  the  synthetic  aperture  L as  caused  by  the  radar 

s 

signal  propagation  through  each  incremental  atmospheric  layer.  This  is 

illustrated  in  Figure  3 for  a layer  dh  thick  at  altitude  h for  a target 

T,  slant  range  R,  and  aircraft  altitude  h . The  effects  of  all  layers 

o 

will  then  be  added  together  to  determine  the  total  atmospheric  effect. 
Porcello  [21  used  this  same  approach.  Strohbein  [7]  used  a vector 
analysis  approach  and  his  results  agreed  with  Porcello's  results. 

Curtis  [8],  in  a related  study,  based  his  analysis  also  upon  a vector 
approach  with  similar  results.  He  then  applied  a separate  model  for 
which  produced  differing  results. 


! 


7.  Strohbein,  J.W.,  Covariance  Functions  and  Spectra  for  Waves  Propa- 
gating in  a Turbulent  Medium  to  or  from  Moving  Vehicles,  IEEE  Trans., 
AP.,  March  1975. 

8.  Curtis,  D.W.,  Atmospheric  Induced  Phase  Error  Effects  in  Synthetic 
Aperture  Radar  Systems,  A EEM  No.  458,  Goodyear  Aerospace  Corporation, 
Litchfield  Park,  Arizona,  July  1973. 


11 


r 


....  — „„ — 


'V.'  * "i. 


Tatarski,  based  upon  his  physical  model  results,  analyzed  the  case 

for  one  way  radio  wave  propagation  along  a path  when  C was  assumed 

n 

constant,  such  as  at  a constant  altitude.  He  examined  the  phase  array 
( structure  function  D^(r)  by  evaluating  the  equation 


Vr) 


R H 2 

= Elk  / n(?,x  + r)  d?  - k I n(?,x)  dC  I 


(5) 


The  geometry  of  the  physical  model  for  this  analysis,  denoted  the  plane 
wave  case,  is  illustrated  in  Figure  4(a).  The  calculated  phase  struc- 
ture function  is  the  expected  value  squared  of  the  propagation  path 
phase  difference  between  two  parallel  ray  paths  over  a total  distance  R 
as  a function  of  path  separation  r.  By  the  use  of  Eqs.  1 and  2 above, 
Tatarski  found  that 


2 2 5/3 

l),(r)  = 8 k R C r 1 

ip  n 


(b> 


where  3 = 2.91. 


The  expression  tor  D^(r)  can  then  be  extended  to  evaluate  the  case 

when  the  value  of  C varies  along  the  path.  That  is  the  case  when  C 
n n 

is  a function  of  the  coordinate  ? in  Figure  4(a).  The  incremental  value 
of  D (n)  due  to  an  elemental  path  length  d?  is 

4> 


2 2 5/3 

d Mr)  = 3 k C (?)  r ' d? 

n 


(7) 


The  total  structure  function  is  therefore 


D (r)  = 3 k‘ 

4> 


K 

/ C" 
r\ 


2 5/3 

(?)  r ' d? 


(8) 


This  result  can  then,  in  turn,  be  extended  to  evaluate  the  SAR 
rase.  The  SAR  case,  which  could  be  called  the  spherical  wave  case,  is 
illustrated  in  Figure  4(b).  The  various  ray  paths  between  a point  being 
imaged  (T)  and  the  synthetic  aperture  (L^)  will  be  examined.  The  desired 
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r 


result  Is  a structure  function  of  the  two  way  phase  along  the  synthetic 

aperture.  The  Incremental  component  of  this  structure  function  with  a 
ri 

ray  separation  of  — is 
K 


d D.(r)  « 4 8 kJ  C2(i) 
ip  n 


and  the  total  structure  function  is 


(*) 


5/3 


d£ 


(9) 


D^(r)  - 4 B k 


(05/3j ■. 

JQ 


2a)  «5/3  de 


(10) 


The  factor  four  Is  due  to  the  two  way  path  factor  that  must  be  included 
for  the  radar  case. 

2 

The  values  of  C are  measured  as  a function  of  altitude  h,  not 
n 

path  distance  Therefore,  a change  of  variables  yields 


V— ^(rfrc> 


,5/3 
h dh 


(ID 


or 


. . 5/3  5/3 

Mr)  * C.  r 


where  C,  A - 

<J>  - h 


4 8k 


h 

f r°c 

ho  J " 


2 (h)  h5/3  dh 
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A SAR  Atmospheric  Phase  Error  Constant 

This  relationship  can  be  further  modified  to  account  for  the  first 
order  effects  of  earth  curvature.  The  resultant  relationship  is 


Vr)  ” 4 p2  k Mo 


2 . [fe(h)  - K3 


5/3 


(j)5/J  °cn2(h)  L 6 . (h  + RJ  dh  (12) 


f (h) 
e 
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„ , . „ 5/3  5/3 

or  Vr)  - s r 


Where  C4>  - R 


2 - f»  2 [fe<">  - K]5/' 

4 S K I c„ <h)  — r® — (" + V dh 

Jo  e 


3/5 


(13) 


K t Tli  2 + 2h  R - R2] 
2 R [_  o o e 


f (h) 
e 


/ 5 2 

vh~  + 2h  R + K 


R = earth  radius, 
e 


The  relationships  and  models  of  Eqs.  12  and  13  are  used  as  the 
basis  of  the  data  analysis  contained  in  the  remainder  of  this  report. 

% 

2.1.3  AVERAGE  SAR  IMPULSE  RESPONSE 

The  average  or  expected  SAR  impulse  response  is  now  determined  as 
a function  of  the  level  of  atmospheric  turbulence.  First,  the  average 
impulse  response  for  normal  fixed  focus  processing  is  derived  assuming 
no  other  sources  of  phase  errors.  Second,  the  response  assuming  a sim- 
ple focus  correction  is  examined.  Finally,  in  this  section,  the  average 
response  is  characterized  by  various  parameters  such  as  measures  of 
mainlobe  width  and  sidelobe  levels. 


2. 1.3. I Fixed  Focus  Average  Response 

The  average  point  target  response  for  a SAR  can  be  determined  as  a 

function  of  the  phase  error  structure  function  D,(r)  as  defined  in  Eqs. 

<P 

11  through  13.  The  received  radar  signals  along  the  synthetic  aperture 
due  to  a single  target  T have  the  fo-  n 


s(x) 


A (x)e 


j<Mx) 


-.1 


2nx 

\R 


(14) 
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The  amplitude  term  A(x)  is  caused  by  antenna  pattern,  radar  range,  and 
receiver  characteristic  effects.  The  random  phase  term  t|>(x>  is  assumed 
to  be  caused  solely  hy  the  atmosphere.  The  remaining  phase  term  is  the 
geometric  variation  in  the  range  to  the  target  along  the  synthetic  array. 
This  is  the  phase  term  which  is  "matched"  in  order  to  create  a radar 
image.  Such  a matched  filter  receiver  can  be  modeled  to  have  an  impulse 
response  of 


h(x)  =*  e 


otherwise 


The  receiver  output  magnitude  is  simply  the  magnitude  of 


A(u)  e e 


The  amplitude  term  over  the"  processing  aperture  can  be  controlled  to 
shape  the  nominal  impulse  response,  especially  to  control  sidelobe 
levels . 

A derivation  of  t he  average  point  target  response  using  results  of 
the  previous  section  is  contained  in  Appendix  A.  The  linear  component 
of  the  phase  error  term  i>(x)  is  removed  while  computing  the  average  re- 
sponse since  that  term  causes  a simple  point  target  image  position  shift 
but  does  not  modify  the  image  shape.  The  resultant  phase  error  is  then 
assumed  to  be  Gaussian.  Results  of  the  derivation  from  Appendix  1 are 
that  the  expected  response  is 


Eifi(x)l2  ■ Ls‘  *1  (jh ; \) 


The  f unc  1 1 on  g 


is  the  normalized  average  point  target  response 
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of  the  radar  and  5 ^ is  the  nominal  unweighted  impulse  response  width  ot 
the  radar  where 


\ IK 


(18) 


This  value  of  3 is  equal  to  half  the  null-to-null  width  of  the  ideal 
unweighted  impulse  response. 

The  parameter  K_  can  he  called  the  synthetic  aperture  phase  error 
coefficient  due  to  atmospheric  anomalies.  The  value  of  this  coefficient 
is  the  product  of  half  the  synthetic  aperture  length  and  the  atmospheric 
phase  error  constant.  It  is  defined  as 


i ^ (!■)> 


N Synthetic  Aperture  Atmospheric 
Phase  Error  Coefficient 


where  C was  previouslv  defined  bv  Eq.  13.  The  value  of  K is  deter- 
mined  by  the  aperture  size,  the  system  geometry,  and  the  atmospheric 
conditions. 

The  expected  variations  of  the  phase  error  across  the  synthetic 
aperture  can  be  determined  from  Eqs.  3,  12,  and  19  to  be 


(C  r) 


5/3, 


(2K  ) 


5/3 


(20) 


r-1. 


r-L 


If  K ”1  then  D is  equal  to  n radians  across  the  aperture.  Table  1 

lists  values  of  D versus  K . The  greater  the  turbulence,  t lie  greater 

<? 

the  value  of  K . When  there  is  no  atmospheric  phase  error  then  K = 0. 

$ ,*> 

The  normalized  average  response  g^s*  ; Kn  j has  been  evaluated. 

The  results  are  illustrated  in  Figure  5 for  five  separate  values  ot  K . 
When  there  is  no  phase  error  and  K » 0 then 
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Figure  5.  Average  Impulse  Resp 
Atmospheric  Phase 
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I'll  In  luilrl  Ion  has  been  normalized  so  tli.it  tin'  peak  value  Is  unity  and 
t lio  I list  mills  ornn  at  x - * 1.  Not  o I rum  Figure  S that  as  tin*  phase 
I'ltof  i'Oi'I  t l c U*nt  Increases  tin*  peak  value  ilivn*tisi*s,  tin*  nulls  begin 
to  till  In,  ami  the  stdolnhe  lev«*ls  l tv  lease. 

The  phase  ei  rot  coettlclent  Is  determined  hv  both  the  synthetic 
aperture  st/e  and  the  atmospheric  phase  erroi  constant.  Thus  the  phase 
error  coot  t l e lent  ean  Inerease  with  either  the  aperture  size  01  atmos 
photic  t urbul  eme . The  eoordluate  x Is  normalized  hv  tin*  apetture  st/e. 
In  other  words,  tor  a tlxed  aperture  the  ettrves  mav  I'*-  compared  dlreet  Iv 
to  II  lust  t ate  the  i-tteet  ol  varying  at  mosphet  l e phase  errors  while  tot 
a varying  aperture  with  t 1 xed  turbulence  the  coordinate  x as  well  as  tin- 
phase  errot  eoel  t I e I out  must  he  varied.  These  general  normalized  curves 
ean  he  used  to  Illustrate  anv  eomhl nation  ol  these  eases  hut  they  must 
he  ear** t ti  I 1 v l n t erpt  ot  »*d  . 

The  I'tiives  ot  Figure  S have  been  redrawn  In  Figure  h to  I I lust  tale 
various  types  ot  ***«mpar  I sons  whteh  eould  he  made.  In  Figure  Mai  t hi* 
normalized  eutves  ol  Figure  have  boon  simply  super  Imposed  t o I I Ins 
(rate  tin*  **t  t eels  ot  varying  l»*v«*ls  ol  turhulenee  with  a tlxed  aperture 
size.  Note  that  as  the  tuihulonee  Ineieases  with  a I l xed  aperture  two 
ma|oi  i-tleets  are  evident.  First,  the  p*>ak  response  deet  eases  and 
seeoud,  the  sldelohe  levels  liu'fi'as*'.  The  total  energy  In  t h«*  t espouse 
st  ays  const  ant  . 

Tin*  eurvi  i are  redrawn  It*  Figure  e*U'^  with  ea*'h  peak  value  set  to 
unity  so  that  the  detailed  struetur**  ean  he  eompared. 

The  curves  In  Figure  h(e1,  |*U,  ami  lei  show  t h«*  data  redt  awn  to 
compute  the  elt**i-t  ot  varying  the  ap«*rtur*‘  sl/e  with  a tlxed  amount  ol 
atmospheric  turbulence.  In  Figure  h|cl  the  i-urves  are  drawn  to  sltv'w 
this  total  etlect.  As  t h»*  aperture  Is  made  larger,  tin*  Impulse  response 
width  narrows,  the  sldelohe  levels  Increase  with  respect  to  the  peak 


.’I 


value,,  and  both  the  total  energy 
curves  are  redrawn  In  Figure  b(d) 
constant  impulse  response  energy, 
for  each  peak  value  set  at  unity. 
1 dB  width,  maximum  peak  value  Is 


and  the  peak  value  increases.  These 
with  the  amplitude  normalized  for 
They  are  also  redrawn  In  Figure  6(e) 
It  can  be  shown  that  the  minimum 
obtained  for  a value  of  K,(>  = 6.7 


These  curves  illustrate  that  the  choice  of  an  optimum  processing 
aperture  must  be  carefully  considered  for  each  application.  Many  times 
the  processing  aperture  is  set  to  yield  a response  close  to  dit fraction 
limited,  l.e.,  with  well  defined  nulls  and  first  sidelobe  levels  close 
to  -12  dB.  This  response  is  generally  thought  to  be  an  optimum  Image. 
The  choice  of  an  optimum  image  under  these  conditions  is  actually  a 
complicated  >ne  and  depends  upon  other  system  parameters,  such  as  Is 
the  syster.i  power  limited  or  signal- to-o lut ter  limited? 


Sidelobe  weighting  can  be  used  to  control  the  impulse  response 

sidelobe  levels  which  are  caused  by  the  processing  aperture  function. 

The  expected  impulse  response  must  be  recalculated  for  such  cases.  I'hls 

2 

has  been  done  for  two  aperture  weightings,  one  a eos“x  and  the  other  a 
30  dB  Taylor  weighting.  The  results  of  these  calculations  are  illus- 
trated in  Figure  7.  These  results  clearly  show  the  effects  of  the 

phase  errors  for  smaller  values  of  K when  the  aperture  fund  Ion  side- 

v 

lobes  are  reduced.  If  a weighting  is  used  to  achieve  low  sidelobes  then 
much  less  phase  error  can  be  tolerated  than  for  the  unweighted  aperture 
case.  Note  that  the  value  of  K,  for  these  cases  was  calculated  tor  6 
as  defined  previously.  In  the  weighted  cases  shown  the  nominal  main- 
lobe  width  is  broadened  by  factors  of  l.H  and  1 . r>  respectively. 


In  summary.  If  the  refractive  index  constant  11  (h)  profile  Is  mea- 
sured then  the  SAR  atmospheric  phase  error  constant  C can  be  calculated 
for  any  system  aircraft  altitude  and  target  range  by  using  Kq.  13.  The 

phase  error  coefficient  K can  then  be  determined  bv  the  processed  aper- 

<P 

tore  size  and  C..  From  the  parameter  K.  and  the  aperture  size  the 
<i)  T 

average  response  is  then  known. 
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(x/d^,  K ) (dB) 


(c)  Varying  Aperture,  Fixed  Turbulence 


Figure  6 


(d)  Varying  Aperture,  Fixed 
Turbulence,  Normalized 
Energy 


Turbulence,  Normalized 
Peak  Values 


Rescaled  Average  Impulse  Response  Plots  (Unweighted) 


2. 1.3.2  Corrected  Focus  Average  Response 

In  the  analyses  summarized  in  Section  2. 1.3.1  and  Appendix  A both 
the  average  value  and  the  linear  component  of  the  phase  error  across 
the  synthetic  aperture  were  ignored  since  they  do  not  affect  the  point 
target  impulse  response  shape.  A third  component  of  the  phase  error  is 
the  quadratic  component  which  causes  a change  in  target  focus.  This 
component  possibly  can  be  removed  by  an  automatic  focusing  processor. 

A derivation  of  the  average  point  target  response  with  the  quadratic 
term  removed  is  contained  in  Appendix  B.  It  is  shown  that 


E| F_(x) | 2 = L 2 g 


^ ^ \) 


2VSo 


(22) 


The  normalized  refocused  response  g„  ( ~ ; K \ has  been  calculated  for 

-\6o  V 

various  values  of  and  is  illustrated  in  Figure  8. 

The  amount  of  focus  correction  or  total  quadratic  phase  error  over 
the  processor  aperture  may  be  found  from  Eq.  B-12  of  the  Appendix  to  be 


/L\2  /L 

s)2 

2 

( ~ 1 O = ( 

i ) 

E(m  ) 

\ 2 / m2  \ . 

! / 

2 

1/2 


= .56  K 


5/6 


(23) 


This  total  quadratic  error  becomes  roughly  1 radian  for  a value  of 


If  the  curves  of  Figures  5 and  8 are  compared  it  is  seen  that  the 
results  are  very  similar.  The  amount  of  focus  correction  is  not  great. 
The  focused  curves  exhibit  a higher  peak  value  and  a narrower  width  main- 
lobe  with  lower  near-in  sidelohe  levels.  The  two  cases,  both  focused 

and  unfocused,  are  illustrated  in  Figure  9 for  the  case  of  K,  = 4. 

<P 

2. 1.3. 3 Impulse  Response  and  Phase  Error  Descriptors 


The  average  impulse  response  pattern  shape  depends  upon  the  level 


of  atmospheric  turbulence,  the  synthetic  aperture  size,  and  the  aircraft 
altitude  and  target  range.  As  shown  in  the  previous  section  the  shape 
is  determined  by  the  phase  error  coefficient  K . It  is  desirable  for 
the  purpose  of  comparing  various  cases  and  for  preparing  system  designs 
to  be  able  to  quantify  the  impulse  response  and  atmospheric  phase  error 
by  simple  parameters  or  analytical  functions.  Five  separate  parameters 
have  been  evaluated  as  possible  descriptors  of  the  impulse  response 


shape.  These  parameters  are  listed  in  Table  2.  They  have  been  calcu- 
lated for  both  normalized  unweighted  responses,  g I K,J  and 

/ \ '° 

®2r^“  ; K^J.  They  are  discussed  below  and  results  are  shown.  Also,  the 
' o ' 

spectral  density  of  the  phase  error  is  presented  for  inclusion  into  a 
phase  error  budget. 


TABLE  2.  IMPULSE  RESPONSE  DESCRIPTORS 


Symbol 

Description 

W 

3 dB  Width 

U 

Equivalent  Rectangle 

I SLR  | 

Sidelobe  Energy 

Total  Energy 

I SLR 

2 

Sidelobe  Energy 

Ma inlobe  Energy 

I SLR, 

3 

Phase  Error  Energy 
Undistorted  Energy 

3 dB  Width 

A common  measure  of  the  resolution  of  a radar  image  is  the  3 dB 
width  of  the  impulse  response.  The  3 dB  width  W Is  defined  bv 

*(*r  ? <0  • 2 8(° : 0 12 

The  value  of  W has  been  calculated  as  a function  of  the  phase  error 
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■ 


— 1 


coefficient  and  is  illustrated  in  Figure  10.  Note  that 
width  increases  slightly  with  K^.  This  measure  of  impulse 
quality  is  nearly  exclusively  determined  by  the  processing 
size,  not  by  the  level  of  phase  error. 


the  3 dB 
response 
aperture 


Equivalent  Rectangle 

The  equivalent  rectangle  p is  defined  as  the  width  of  a rectangle 
whose  height  is  equal  to  the  average  response  peak  value  and  whose 
energy  is  the  same  as  the  average  response  energy.  This  definition  can 
be  expressed  as 


P 


dx 


(25) 


which  becomes 


P 


RX 

2L 

s 


(26) 


If  there  are  no  atmospheric  phase  errors  then  P equals  . The  value 

Js 

of  P is  simply  a measure  of  the  peak  value  of  the  average  response  and 
therefore  is  not  a very  useful  direct  measure  of  image  quality.  The 

dependence  of  p on  K is  shown  in  Figure  10.  For  the  case  analyzed  in 

9 

this  study  the  decrease  in  peak  value,  or  increase  in  p,  is  most 
directly  related  to  the  increase  in  impulse  response  sidelobe  levels. 
Therefore,  the  parameter  p is  more  nearly  a measure  of  sidelobe  levels 
rather  than  mainlobe  width. 


ISLRi 

One  measure  of  integrated  sidelobe  energy,  ISLRj , is  defined  as 
the  sidelobe  energy  of  the  normalized  impulse  response  divided  by  the 
total  response  energy.  For  this  purpose  the  separation  point  between 
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the  malnlobe  and  sidelobe  is  defined  as  the  point  where  the  first  zero 
would  occur  in  the  pattern  if  there  were  no  phase  errors.  ISLR^  can  be 


written 


-i 

li  = j s(x^ 


) dx  + / g(x;K^)  dx 


1 

- J g(x; 


K.)  dx 

4> 


The  dependence  of  ISLR^  on  K is  also  shown  in  Figure  10. 


A closely  related  measure  of  sidelobe  energy,  ISLR^,  is  the  ratio 


of  sidelobe  to  mainlobe  energy.  It  can  be  written  as 


y g(x;K^)  d x + y g(x;Kj,)  dx 


islr2  = 


g(x;K  ) dx 


1 - ISLR, 


This  is  also  illustrated  in  Figure  10. 


A third  measure  of  sidelobe  energy,  ISLRy  is  the  energy  contained 
in  the  impulse  response  due  to  the  distorted  phase  divided  by  the  energy 
due  to  the  undistorted  phase.  This  definition  was  used  by  Curtis. 

When  the  phase  errors  are  small  we  can  assume  that 


r 
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eW(X>  ; I + W(X) 


The  impulse  response  from  Eq . 16  is 


F (x)  = 


/%„ 


) [ 1 + j4>  (u)  ] e 


The  total  response  energy  for  A(u)  = 1 is 


W 


2 R\  / L 9 

F(x)|  dx  = — 1+1  |<Ku)!  “ du 


Therefore 


islr3  = r 


f/ 


i 4> (w)  | du 


This  has  been  evaluated  and  is  also  shown  in  Figure  10.  Table  3 listing 
values  for  ISLR  has  been  prepared  since  the  curves  cannot  be  easily 
read  for  small  values  of  K.. 


Spectral  Density 


The  spectral  density  of  the  total  phase  error  function  is  an  impor- 
tant parameter  to  a system  designer.  Equations  3 and  13  showed  that 

D (r)  - C 5/3  r5/3  (33) 

<p  <p 

= E[ip(x  + r)  - $(x)l2 
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It  can  be  8 hewn  that  the  spectral  density  of  this  phase  error  along  the 
aperture  Is 


2.2  EXPERIMENTAL  DATA 

Experimental  measurements  of  the  atmospheric  Index  oi  refraction 

were  made  by  the  U.S.  Department  of  Commerce.  The  data  collected  wore 

reduced  into  profiles  of  the  refractive  index  constant  C (h).  These  ro- 
il 

suits  were  then  provided  to  ERIM  for  analysis  of  t lie  expected  SAR  per- 
formance as  a function  of  geometric  radar  parameters  based  upon  the 
model  described  In  Section  2.1.  The  results  of  this  data  collection 
and  analysts  are  described  in  this  section.  In  Section  2.2.1  the  data 
collection  phase  is  briefly  reviewed.  In  Section  2.2.2  the  data  analy- 
sis results  are  presented. 


2.2.1  DATA  COLLECTION 

A separate  experimental  aircraft  flight  test  program  was  carried 

out  by  the  U.S.  Department  of  Commerce,  Office  of  Telecommunic.it  ions . 

The  purpose  of  this  flight  test  program  was  to  make  measurements  of 

index  of  refraction  profiles  to  be  used  in  this  analysis  of  the  effect 

of  atmospheric  variations  upon  synthetic  aperture  radar  performance. 

The  data  collected  was  reduced  into  profiles  of  the  index  ot  refraction 
2 

constant  C^  and  was  then  supplied  to  ERIM  for  t lie  appl  icat  ion  to  SAR 
performance  estimates.  The  data  collected  and  supplied  to  ERIM  are 
presented  in  detail  in  a separate  report  |6]  published  bv  the  D.S. 
Department  of  Commerce  and  are  only  briefly  presented  herein. 

The  data  collection  was  based  upon  measurement  techniques  and  equip 
ment  which  have  been  developed  over  a period  of  several  years  at  the 
U.S.  Department  of  Commerce  [R],  The  instrumentation  package  is 


‘1 . Vetter,  M..I.  anil  M.C.  Thompson,  Jr.,  Solid-State  Microwave  Refraeto- 
meter,  Rev.  Sci.  Instrum.,  Vol.  38,  No.  12,  1967,  pp.  1726-1727. 
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relatively  small  (total  weight  11  pounds)  and  can  be  rapidly  Installed 
and  used  on  a wide  variety  of  small  aircraft.  The  major  component  of 
the  equipment  Is  an  ITS  model  h X-band  ref ractometer  which  provides 
analog  measurements  of  the  atmospheric  Index  of  refraction,  temperature, 
and  pressure.  In  the  current  system  these  measurements  are  filtered, 
digitized,  and  recorded  on  magnetic  tape  for  subsequent  analysis.  The 
data  collection  system  was  modified  specifically  for  this  program  to 
provide  an  Increased  spatial  data  sampling  rate. 

During  i his  program,  three  sets  of  data  gathering  flights  were 
flown  for  a total  of  17  separate  flights.  The  Initial  set  ol  five 
flights  were  flown  In  an  area  near  t hi'  boulder,  Colorado,  11. S.  Depart- 
ment ol  Commerce  laboratory.  These  Initial  flights  were  carried  out 
for  the  purpose  of  testing  and  calibrating  the  modified  Instrumentation 
package  and  to  provide  a sample  data  base.  The  second  set  of  nine 
flights  were  flown  In  an  area  close  to  I'.gltn  Mr  Force  base,  Florida. 

The  purpose  of  those  flights  wire  to  gather  data,  for  the  analysis, 
which  would  exhibit  samples  from  relatively  severe  atmospheric  conditions. 
The  Eglin  area  was  chosen  due  to  Its  close  proximity  to  a land/water 
boundary  and  due  to  the  prevailing  local  weather  (temperature/humldltv) 
conditions.  A third  set  of  three  additional  Mights  were  flown  near 
the  Boulder  laboratory.  The  purpose  of  these  final  flights  was  to  pro- 
vide additional  data  up  to  higher  altitudes  during  a different  season 
of  the  year.  The  total  amount  of  data  gathered  were  severely  limited 
in  terms  of  numbers  of  flights,  seasonal,  diurnal,  and  geographic  con- 
siderations but  were  deemed  sufficient  to  provide  significant  results 
from  the  subsequent  analysis. 

Two  different  types  of  flight  paths  were  flown.  One  type  used  was 
composed  of  a series  of  level  runs  at  1,000  ft  altitude  intervals.  The 
other  type  was  a steady  and  shallow,  spiral  or  straight,  ascent  or 
descent.  The  specific  flights  carried  out  are  summarized  In  Table 
Flights  were  flown  to  altitudes  as  high  as  29,000  It.  Data  were 
gathered  under  a variety  of  weather  conditions  although  all  flight 
patterns  were  chosen  so  that  VFK  conditions  were  maintained.  Data  were 

IS 


TABLE  4.  FLIGHT  SUMMARY 


TLT. 

WO. 

0 ATT /TIME 

LOCATION 

MAX. ALT. 

x-rt. 

TYPE  REMARKS 

1 

8/30/76 

1000/1130 

Colorado 

9.5 

* Ladder  ascent  with  level  runsj  constant 
descent  on  E,  W course  from  Longmont 

2 

8/3/76 

1300/1500 

- 

12.5 

mm  m m m m m 

3 

9/21/76 

0830/1130 

r. 

12.5 

mm  m m m m m 

4 

9/28/76 

1230/1410 

« 

11.5 

Slant  aacont  from  Longmont  towards  Jackson 
Raaervolr;  ladder  descent  with  level  runs 
on  NE,  SW  course  from  Jackson  Reservoir 

S 

9/29/76 

1230/1530 

■ 

12.5 

« m n n n n n 

6 

10/8/76 

1600/1830 

Florida 

20.5 

Spiral  ascent,  ladder  descent  from  Panama 
City  towards  Dothan,  ALA. 

7 

10/11/76 

0930/1230 

■ 

23.0 

Spiral  ascent,  ladder  descent  takeoff 
Dothan,  ALA,  spiral  up  east  of  airport. 

• 

10/12/76 

0840/1050 

• 

23.3 

Spiral  ascent,  ladder  descent  over  coast- 
line - 10  mi.  S of  Tyndall  AFB 

9 

10/13/76 

2100/2400 

■ 

23.6 

Slant  ascent  east  and  descent  west, 

Panama  city  to  Jacksonville 

10 

10/14/76 

1330/1630 

■ 

24.0 

Slant  ascent  and  descent  over  Gulf  on 

190*  radial  from  Panama  City 

11 

10/15/76 

1000/1400 

• 

25.0 

Slant  ascent  east  to  Jacksonville,  descent 
west  to  Panama  City 

12 

10/16/76 

1030/1340 

« 

25.0 

Slant  ascent  SE  towards  Tampa,  descent 

MM  to  Panama  City,  FLA,  following  the 
coastline  on  both  ascent  and  descent - 

13 

10/18/76 

1100/1400 

m 

25.0 

Slant  ascent  and  descent  over  Gulf  on 

185*  radial  from  Panama  City,  FLA 

14 

10/19/76 

1000/1300 

m 

24.0 

Slant  ascent  toward  west  to  Gulfport,  MISS, 
descent  east  to  Panama  City,  FLA,  following 
coertline  on  both  ascent  and  descent. 

IS 

11/16/76 

1430/1700 

Colorado 

28.5 

Ascent  and  descent  “figure  8"  vertical, 
end  Longmont,  CO,  as  turning  points 

16 

11/17/76 

1230/1530 

m 

28.0 

Slant  ascent  toward  north  from  Boulder, 
descent  south  toward  Boulder 

17 

11/19/76 

0900/1230 

■ 

29.0 

Slant  ascent  toward  east  from  Boulder. 

descent  west  toward  Boulder 
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gathered  during  rain  and  on  a few  occasions  raindrops  caused  abnormal 
refract ivltv  readings  so  that  data  collected  at  these  times  were  ignored. 

The  measured  Index  of  refraction  (n)  data  were  used  to  calculate 

2 

index  ot  refraction  structure  function  (C  ")  profiles.  Figure  II 

illustrates  examples  of  typical  measured  index  of  refraction  data.  The 

data  were  measured  such  that  independent  samples  were  obtained  at  a 

2 

spatial  rate  of  roughly  one  per  meter.  F.ach  of  the  C profiles,  along 

n 

with  temperature  and  refractivity  profiles  are  presented  in  the  U.S. 

<7  2 

Department  ot  Commerce  report.  The  C profiles  are  reproduced  in 

2 n 

Figure  12.  The  values  of  C^~(h)  were  computed  at  1,000  ft  altitude 
intervals  for  each  data  set  except  sets  5-D,  6-D,  and  7-D  where  500  ft 
Intervals  were  used. 

The  majority  of  the  atmospheric  water  vapor  which  causes  the  sub- 
ject index  of  refraction  variations  Is  distributed  in  the  lower  15,000 

ft  of  tlie  atmosphere.  The  measured  data  showed  that  principal  values 
2 

of  C were  below  the  10,000  ft  altitude.  In  a few  cases,  there  were 
n 

significant  values  up  to  20,000  ft  but  in  all  of  these  cases  this 
amounted  to  a fraction  of  the  peak  value. 

2.2.2  DATA  ANALYSIS 

This  section  contains  the  results  of  the  data  analysis  which  was 

2 

carried  out  at  ERIM.  First,  a brief  analysis  of  the  C profiles  which 
was  carried  out  to  assess  the  validity  of  the  analytical  atmospheric 
model  is  presented.  Then  the  results  of  an  in-do.pth  data  analysis  to 
predict  SAK  performance  for  most  of  the  test  flights  is  described. 
Finally,  these  results  are  summarized. 

2.2. 2.1  Analytical  Model  Examination 

Refraction  index  constant  profiles  were  calculated  at  the  U.S. 
Department  of  Commerce  from  the  index  of  refraction  measurements.  The 
calculations  performed  are  derived  from  Eqs.  1 and  2 to  bo 
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V 

FLIGHT  3,  DESCENT  (NO  DATA) 


Figure  12.  Refractive  Index  Constant  Profiles  by  Flight  Number 
and  Ascent  (a)  or  Descent  (d)  (continued) 
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Figure  12.  Refractive  Index  Constant  Profiles  by  Flight  Number 
and  Ascent  (a)  or  Descent  (d)  (continued) 
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Figure  12.  Refractive  Index  Constant  Profiles  by  Flight  Number 
and  Ascent  (a)  or  Descent  (d)  (continued) 
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Refractive  Index  Constant  Profiles  by  Flight  Number 
and  Ascent  (a)  or  Descent  (d)  (continued) 
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Figure'  12.  Refractive  Index  Constant  Profiles  by  Flight  Number 
and  Ascent  (a)  or  Descent  (d)  (continued) 
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Figure  12.  Refractive  Index  Constant  Prof l lea  by  Flight  Number 
and  Ascent  (a)  or  Descent  (d)  (continued) 
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Refractive  Index  Constant  Profiles  by  Flight  Number 
and  Ascent  (a)  or  Descent  (d)  (continued) 
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Figure  12. 


Refractive  Index  Constant  Profiles  by  Flight  Number 
and  Ascent  (a)  or  Descent  (d)  (continued) 
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Figure  12.  Refractive  Index  Constant  Profiles  by  Flight  Number 
and  Ascent  (a)  or  Descent  (d)  (continued) 
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(35) 


C 2(h)  = E[n(x  + r)  - n(x)]2 

r ~ 


These  profiles  were  calculated  for  1,000  ft,  or  on  one  occasion  500  ft, 

altitude  intervals.  As  a test  of  the  atmospheric  turbulence  model  the 

2 

value  of  was  computed  for  values  of  path  length  r ranging  between 
10  and  500  meters.  The  calculated  values  must  be  independent  of  r and 
dependent  only  upon  h in  order  for  the  model  to  be  valid.  Figure  13 
shows  two  samples  of  these  calculations.  It  is  seen  that  the  dependence 
upon  the  path  length  is  indeed  minimal  and  accordingly  the  model  seems 
to  be  validated. 

2. 2. 2. 2 SAR  Performance  Analysis 

It  was  shown  in  Section  2.1.3  that  the  normalized  impulse  response 
shape  is  characterized  by  the  synthetic  aperture  phase  error  coefficient 
K. . From  Eqs.  18  and  19  this  coefficient  may  be  written  as 


where  = synthetic  aperture  phase  error  coefficient 

K,k  = normalized  coefficient  for  (S  = 1 ft 
H 1 o 

C = SAR  atmospheric  phase  error  constant 
5 = nominal  unweighted  impulse  response  width 

X = radar  wavelengtii  (assumed  0.1  ft) 

R = slant  range 

2 

The  measured  C (h)  profiles  were  used  to  calculate  the  constant 
n 

C\, . The  relationship  of  Eq.  13  was  applied  for  various  values  of  slant 
range  R and  altitude  h.  Then  Eq.  36  was  used  to  calculate  values  of 
. These  values  were  calculated  for  ranges  out  to  100  nmi  and  for 

altitudes  of  50,  70,  and  90  thousand  ft.  Results  of  this  data  analysis 
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are  presented  In  Figure  14.  The  value  of  wavelength  used  was  0.1  ft 
and  the  4/3  earth  radius  model  was  used. 

These  curves  were  calculated  for  each  data  set  for  which  complete 
data  were  available.  Linear  interpolation  was  used  to  replace  data 
values  missing  due  to  radio  interference  problems.  The  effect  of  this 
replacement  should  be  negligible.  No  attempt  was  made  to  use  data  with 
values  missing  due  to  rain  saturation  of  the  measuring  instruments. 

These  curves  can  be  used  to  directly  estimate  SAR  imagery  perfor- 
mance. The  image  coefficient  Kt  can  be  directly  obtained  for  any  flight 
and  any  set  of  system  parameters;  range,  altitude,  and  nominal  impulse 
response  width.  As  an  example,  suppose  the  system  design  is  R = 50 
miles,  h = 50,000  ft,  6 = 10  ft,  then  for  flight  2-A  the  coefficient 
K„  is  simply  approximately  ~ x 1.1  or  0.11.  Thus  the  effect  would  be 
negligible  even  for  the  weighted  case.  For  the  worst  case,  flight  1 ID 
the  same  case  yields  = 0.42.  For  the  same  flight,  if  the  range  was 
extended  to  100  nmi  the  K,  = 1.6. 

4> 

In  an  attempt  to  reduce  the  data  even  further  three  performance 
frequency  plots  were  generated  from  the  experimental  data  as  shown  in 
Figure  15.  Each  of  the  curves  represents  a quartile  of  the  data  com- 
puted from  the  14  useable  test  flights.  It  should  be  emphasized  that 
these  plots  are  derived  from  limited  data  and  cannot,  in  general,  be 
used  to  predict  specific  SAR  performance. 

Curtis  [8]  uses  a model  that  predicts  SAR  performance  when  surface 
water  vapor  pressure  is  known.  The  data  analyzed  herein  was  also 
examined  to  assess  that  model.  Results  were  taken  for  an  altitude  of 
50,000  ft  and  a range  of  100  nmi.  The  scatter  plot  obtained  is  shown 
in  Figure  16,  each  of  the  14  Florida  data  flights  is  represented  on 
this  plot.  The  two  quantities  appear  to  be  unrelated.  The  correlation 
coefficient  was  calculated  to  be  0.25. 

From  Eq.  12  it  was  shown  that  the  factor  effecting  performance  was 
2 5/3 

C (h)  h . Surface  measurements  probably  should  not  be  used  alone. 
n 2 

As  C (h)  varies  with  altitude,  the  higher  altitude  values  receive  a 
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COLORADO 
FLIGHT  2-A 


FLIGHT  2,  DESCENT  (NO  DATA) 


Figure  14.  Normalized  Synthetic  Aperture  Phase  Error  Coefficient, 

, for  Experimental  Data 
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FLIGHT  4,  ASCENT  (NO  DATA) 
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FLIGHT  4-D 


Range  (nml) 


Figure  14.  Normalized  Synthetic  Aperture  Phase  Error  Coefficient 
K,  , for  Experimental  Data  (continued) 
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Figure  14.  Normalized  Synthetic  Aperture  Phase  Error  Coefficient 
K.  , for  Experimental  Data  (continued) 
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heavy  weighting.  The  experimental  data  do  exhibit  significant 
disturbances  at  altitudes  up  to  10,000  ft. 
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Figure  15.  Frequency  Plots  for  6 Vs.  Range,  Altitude 
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CONCLUSIONS  AND  KKCOMMKNDAT 1 ONS 


Wlii lo  the  quality  ol  the  data  available  for  this  study  is  excellent, 

2 

tilt'  quantity  is  limited.  More  T (It)  data  are  needed;  these  should  he 

it 

collected  in  both  Florida  and  Colorado,  as  well  as  in  many  other  geogra- 
phic. areas.  Additional  data  from  both  Florida  and  Colorado  could  he 

•1 

used  to  determine  both  the  spatial  and  temporal  variability  ol  C ' (h). 

The  average  resolution  computed  t rom  the  Colorado  data  was  I liter  than 
from  the  Florida  data,  indicating  that  geographic  factors  should  he 
considered  when  predicting  resolution.  There  are  not  enough  data  avail- 
able to  determine  which  are  the  important  factors  and  how  they  relate 
to  the  limiting  resolution.  The  ultimate  goal  here  is  to  achieve  a 
model  for  the  atmosphere  which  applies  for  many  areas  ot  the  earth. 

An  effort  is  needed  to  relate  computed  resolution  to  available 
weather  data.  It  does  not  appear  that  resolution  can  he  predicted  re- 
liably when  only  surface  weather  data  are  used.  It  may  be  that 
resolution  is  more  closely  related  to  temperature  and  pressure  pro 

files.  These  profiles  were  measured  as  the  C ‘(h)  data  were  gathered, 

n 

although  their  use  as  predictors  was  not  attempted.  it  appears  that 
smooth  profiles  and  good  resolution  occur  together,  but  that  a ragged 
profile  does  not  necessarily  imply  poor  resolution. 

The  atmospheric  model  used  for  this  study  was  Tatarski's  turbulent 
medium  theory.  This  theory  should  be  extended  It'  include  the  effects 
ot  tain.  in  addition,  the  effect  ot  large  amplitude  changes  in  the  data 
should  be  studied.  Those  changes  are  much  larger  than  anv  changes  pre- 
• I . t e.l  bv  t he  mode  1 . 

. . . wet  e no  data  gathered  within  heavy  cloud  formations,  and  it 
•wn  whet  bet  this  model  won  hi  apply  in  such  cases.  Data  within 
, .elicit  .si  and  the  usefulness  ot  the  model  should  be 


APPENDIX  A 

AVERAGE  POINT  TARGET  RESPONSE  FIXED  FOCUS 

Tills  appendix  contains  the  derivation  of  the  average  point  target 
a.-  limit  h response  pattern  tor  a synthetic  aperture  radar  In  the  presence 
ot  atmospheric  phase  errors. 

The  radar  response  to  a point  target  In  the  azimuth  coordinate  is 
modeled,  from  Eq.  16,  as 


l. 


f(x)-  / 


, A ii  xu 

(u)  J XR  . 
A(u)  e e du 


(A—  1 ) 


The  linear  component  of  ^(u)  does  not  affect  the  shape  of  the  response, 
only  Us  location.  Therefore,  the  phase  term  will  be  expanded  by 

>1(11)  “ <Mu)  - m u (A-2) 

and  the  shape  of  the  pattern  will  then  be  determined  by  evaluating 
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<J>(u)  du  + / du  (A-4) 
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Since  the  above  sum  is  non-negative  then  the  minimum  occurs  when 
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12 
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s -L 


4>(u)  du 
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The  factor 
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4>  ( 0)  u du  = 0 


s -L 


can  be  added  to  A-5  to  yield 
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[ 4>  (u)  - ( 0)  ] du 


The  average  point  target  response  is  then 
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(A-5) 


(A-6) 


We  now  use  the  assumption  that  Y^(u)  - Y^(v)  is  Gaussian  to  write 


jlY1(u)-Y1(v)] 


hlY.tuKWl2  -!»v  (u-v) 

1 1 1 = e * '1 


(A-8) 


so  that 
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X(u-v)  1 (Lsu  LsW\ 


du  dv 


(A-9) 


Next  to  find  the  structure  function  of  Y-^(x)  in  terms  of  the  structure 
function  of  4>(x)  we  write 


I)  (x-y)  = E 


= E 


hj (x)  - Y1(y)l‘ 


I4>(x)  - *(y)]‘ 


- 2 (x-y)  E 


n^IiMx)  - 4> (y ) 1 


+ (x-y)  E 


m. 


(A- 10) 


These  three  terms  of  A-10  will  be  evaluated  separately. 

.21 


First,  E 


[if(x)  - <t>(y) I 


D.  (x-y) 


(A- 11) 
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Second,  the  identity 


(a-b)  (c-d)  = |[(a-d)2  + (b-c)2  - (a-c)2  - (b-d)  ] 


(A-12) 


will  be  used  to  write 
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Finally, 


(A-13) 


continued  on  next  page 
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J J xy  EjO(x)  - 0(0)]“  + [0(o)  - 0(y)j 


s L L 
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— [0 ( x ) - 0(0)]  | dx  dy 


ff-fy-h) 
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When  the  three  terms  of  A-10  are  combined  the  result  is 


I\.  (x-y)  = D^(x-y) 


(x-y) 


^ "["if  ~y)  - \{f  -») 
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du  dv 


s <1  •'-l 


This  expression  can  be  simplified  by  using  the  model  for  D^(,r) 


Vrt-c/'3  w5/3 

0 0 


Then  Eq.  A- 15  becomes 
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The  two  integrals  above  can  be  evaluated.  The  first  is 


(A-17) 
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P 

(A-18) 


The  second  is 


.764 


Thus  A-17  becomes 
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where  q(u,v) 
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The  average  point  target  response  can  thus  be  written,  from  Eqs. 
A- 9 and  A-20,  as 


E|f(x)|2  ' Ls2  g(^f-  ; 


(A-22) 


where  g(x;  ) - 


e *"  ^ du  dv 


(A- 23) 


If  we  set  the  weighting  A(u)  = 1 then 


E|  f(x)|2  = Lg2  g 


f!V . A 

1 \ Ar  ’ K4>  / 


(A-24) 


where  g^(x; 


1 f1  r1  -jTT(u-v) 
-1  -1 
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(A-25) 


The  function  g, (x;  K ) has  the  following  properties: 
1 


, / sin  ttx  \ 4 

8l<x;  °> 


(/1-26) 


and 


fn*  s 


) dx  = 1 


(A-27 ) 


The  function  has  been  evaluated  for  a range  of  values  of  x between  0 
and  6 and  for  values  of  K between  0 and  6.7.  The  results  are  shown  in 
Figures  5 and  6. 
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APPENDIX  B 

AVERAGE  POINT  TARGET  RESPONSE  CORRECTED  FOCUS 


This  appendix  contains  the  derivation  of  the  average  point  target 
azimuth  response  pattern  for  a synthetic  aperture  radar  in  the  presence 
of  atmospheric  phase  errors  with  optimum  focussing. 

The  analysis  in  this  appendix  is  very  similar  to  that  in  the  pre- 
vious appendix.  The  phase  term  to  be  expanded  is 


Y2(u)  = <(>  (u) 


m^u  - m^u 


(B-l) 


To  find  m^ 


and  m 


2 


to  minimize 


(B-2) 
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The  minimum  value  is  when 


m,  = — ^ I ucp(u)  du 
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and 


so  [2  : 
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4>(u)  du 


As  before  we  will  evaluate  D (x-y). 
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These  three  terms  must  be  evaluated.  The  first  was  previously  evaluated 
in  Appendix  A.  The  first  portion  of  the  second  term  is 
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(B-6) 


since 


E[  4>(x)  - 4>(y)  ] = 0 


(B-7) 


Using  the  identity  in  Appendix  A,  this  becomes 
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The  remainder  of  the  second  term  is 


Ej2(x2-y2)  m2m1(x-y)^  = 0 


(B-10) 


since  Elm.^]  = 0 


For  the  third  term 
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Since  E[0(x)  - 0(0)]  = E[0(y)  - 0(0)1  = 0 


The  algebraic  identity  yields 


+ [0(0)  - 0(y) ] 2 - [ 0 (x ) - 0 (y) ] 2 | dx  dy 

(continued  on  next  page) 
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After  combining  all  terms  we  have 
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D (x-y)  = D (x-y)  - 2(x  -y  ) 
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By  use  of  D (r)  - C 5/3  | r | 5/ 3 we  find 
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Both  integrals  can  be  evaluated  and  the  average  pattern  with  focus 
correction  is 


E | f (x) | 2 = Lg2 
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where  h(u,v)  = q(u,v)  - - (u2-v2)  j^| u | 5 ^ 3 - | v | 5 ^ 3 J 
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To  write  E | f (.x)  | when  A - 1 in  normalized  form  requires 
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(B-18) 
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Figures  8 and  9 show  g2(x;  K^)  for  x between  0 and  6 and  varying  values 
of  y. 
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